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Abstract 

In this work we reinvestigate braneworld models with nonstandard dynamics. We choose a 
specific model diffused in the literature and focus our attention on the matter energy density, the 
energy of system, the Ricci scalar and the thin brane limit. As the model is classically stable and 
capable of localize gravity, as a natural extension we address the issue of fermion localization of 
fermions on a thick brane constructed out from one scalar field with nonstandard kinetic terms 
coupled with gravity. The contribution of the nonstandard kinetic terms in the problem of fermion 
localization is analyzed. It is found that the simplest Yukawa coupling ij^/^ support the localiza- 
tion of fermions on the thick brane. It is shown that the zero mode for left-handed can be localized 
on the thick brane depending on the values for the coupling constant r/. Furthermore, we present 
an argument that could jeopardize the results of previous work on massive modes. 
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I. INTRODUCTION 



In the last decade, the braneworld scenario has attracted a lot of interests for it gives an 
effective way to solve the hierarchy problem by introducing two 3-branes which are embedded 
in a five-dimensional anti-de Sitter (AdS^) space-time As another attractive property, 
the Newtonian law of gravity with a correction is also given in this braneworld scenario J2]. 

n n 

In the Randall-Sundrum model [l], we can further add scalar fields 3| with usual dynamics 
and allow them to interact with gravity in the standard way. In this scenario, the smooth 
character of the solutions generate thick brane with a diversity of structures In the 

braneworld scenarios, an important issue is how gravity and different observable matter 
fields of the Standard Model of particle physics are localized on the brane. It has been 
shown that, in the Randall-Sundrum model in 5-dimensional space-time, graviton and spin 
field can be localized on a brane with positive tension j^J, sj-joj. Moreover, spin 1/2 and 
3/2 can be localized on a negative-tension brane [9fl. The localization problem of spin-1/2 

nn 

fermions on thick branes is interesting and important [8J-|30[. In order to achieve localization 
of fermions on a brane with positive tension, it seems that additional interactions except 
the gravitational interaction must be including in the bulk. 

On the other hand, the first recent observations 3JJ- 32J have led us with the intriguing 



fact that the Universe is presently undertaking accelerated expansion. These information 
directly contributed to establish some important advances in cosmology, one of them being 
the presence of dark energy. The presence of dark energy has opened some distinct routes of 
investigations. In recent years, there appeared some interesting models with noncanonical 



dynamics with focus on early time inflation or dark energy [33lj - 36], as for instance, the 
so-called k-fields, first introduced in the context of inflation |36j and the k-essence models, 
suggested to solve the cosmic coincidence problem 35(, 37]]. The interaction between dark 



energy and fermion fields has already a precedent in the cosmology context [38]. In the 
context of braneworld scenario, the effect of general brane kinetic terms for bulks scalars, 
fermions and gauge bosons in theories with and without supersymmetry has already been 



analyzed in 



39] . We believe that the conditions for obtaining normalizable zero modes in a 



brane model with generalized dynamics deserves to be more explored. 

In this paper, we reinvestigate a braneworld model with nonstandard dynamics. The 
model C = K{X) - V (</)), where K (X) = X + a\X\X (type I model in Ref. [33]) is 
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considered. We will focus our attention mainly on the matter energy, the energy of system, 
the Ricci scalar and the thin brane limit. As the model is classically stable and capable of 
localize gravity, additionally we address ourselves to the issue of fermion localization on a 
thick brane constructed out from one scalar field with nonstandard kinetic terms coupled 
to gravity. We use the analytical expressions for small a and investigate the contribution 
of this nonstandard kinetic terms in the problem of fermion localization. We find that the 
simplest Yukawa coupling r]^/<j)^/, where rj is the coupling constant, which allows left-handed 
to posses a localized zero mode on the thick brane under some conditions on the value for the 
coupling constant rj. For a, which is not necessarily small, we can not get any expressions 
for the solution of this model, for this case the numerical study is used. The numerical 
results bear out our results for small a. The organization of this paper is as follows: in Sec. 
II, we review and finish the analysis of brane model with generalized dynamics developed 
by Bazeia and collaborators 33j. In Sec. Ill, we study the localization of spin-1/2 fermion 
for this model, we analyze the essential conditions for the localization with the simplest 
Yukawa coupling and present an argument that could jeopardize the results of previous 
work on massive modes. Finally, our conclusions are presented in Sec. IV. 



II. REVIEW OF SYSTEMS WITH GENERALIZED DYNAMICS 



The action for this kind of system is described by 33 ] 



S 



d 5 x\/\g\ 



--R + C{<f>,X) 



(1) 



where g = Det(g a b) and X — | V a 0V a 0. The line element of the five-dimensional space-time 
can be written as 

ds 2 = g ab dx a dx b = e 2A{y) Vflu dx^dx u - dy 2 , (2) 

where we are using the five-dimensional Newton constat AnG^ = 1, y = x 4 is the extra 
dimension (the Latin indices run from to 4), rj^ is the Minkowski metric with signature 
(+, — , — , — ) and e 2A is the so-called warp factor (the Greek indices run from to 3). We 
suppose that A = A{y) and = <fi(y). 

One can determine the static equations of motion for the above system are of the form 



(C x + 2XC XX ) 0" - (2XC X(P - £ ) = -AC X A'. 



(3) 
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A" + 2A l2 = \c, (4) 

A' 2 = l -{C-2XC X ) 1 (5) 

where prime stands for derivate with respect to y, Cx = dC/dX and = dC/dcj). 
Furthermore, the matter energy density is given by 

p(y) = -e 2A ^C, (6) 

and the scalar curvature (or Ricci scalar) is given by 

R = -A(hA' 2 + 2A"). (7) 

The Lagrangian density C(<p, X) has the form 

C = K(X)-V(cj>), (8) 

where K(X) and V(4>) are nonstandard kinetic term and potential respectively. 

For this case, from Eqs. (J3J), (@| and (jSJ) the equations of motion can be expressed as 

(K 1 + 2XK") (j)" -V+ = -AK'A'fi, (9) 

A" + 2A' 2 = -(K-V), (10) 
3 

A' 2 = - (K-V-2XK'). (11) 
3 

These equations are the static equations of motion of a system with nonstandard dynamics. 
In 33), the authors present two explicit models for K(X), here we review one of these 
models. 



A. The model: K(X) = X + a\X\X 



This model is also considered in 40], where a is a real non- negative parameter and 



X = — \4>' 2 - If a = the standard scenario is restored. For this model, the equations of 
motion are 

0" + AA'fi -V <t> = -a (30" + Acj)' A') <P' 2 , (12) 
A" + 2A' 2 = - 1 - (l + | <P' 2 )<p> 2 - 2 -V, (13) 
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A ' 2 = l( 1 + l a( f ) ' 2 )^ 2 -l v - ( 14 ) 



It is possible to rewrite ( 113]) and ( 114]) as 

A" = -^0' 2 (l + a0' 2 ), (15) 
Now, to extend the first-order framework to the braneworld scenario, we follow the work 



in Ref. and choose the derivative of the warp factor with respect to the extra dimension 
to be a function of the scalar field 



A' = - l -W{4>). (16) 



Substituting f[T6]) into f|T5l) . we get 



0' + «0 /3 = i^, (17) 

this equation is a cubic equation in 0', then the real solution to this cubic equation is given 
by 

^=—a (18) 

where 

1 /3 

m(W$) = [BWWj, + 6^ (16a 3 + 27a 4 W$) 1/2 ) . (19) 

It is instructive to note that the equation f|T8|) is the first-order differential equation for the 
scalar field <fi. 

The potential is obtained by substituting ( ITBl) in f fT4j) 



K(0) = i0' 2 + ^0' 4 -^(0) 2 , (20) 



where <f)' is given by f fT8|) . 

On the other hand, we consider the energy functional 4\ 



E[A,(f>] = j dy(-£ system ), (21) 
where C sys t e m = \f\g~\ [— -R/4 + £(</>, X)]. This energy functional is associated with the 



energy of system, as was reported by Hawking [41]. In our case, (|2T|) becomes 



E[A.o} = j ( U,^{\<i>' 2 -ZA^ + ^^ + V\. (22) 



The Euler-Lagrange differential equations from the functional (1221) are the static equations 
of motion f JT2|) and (TIB"]) , this result is valid for all a. Otherwise, substituting the potential 
([20]) in (I22p . for a which is not necessarily small, we can not deduce the first order differential 
equations ([TBI) and f[T8"l) from the energy functional, therefore the solutions of the first order 
differential equations not necessarily minimize the energy functional. 

At this point, it is also instructive to analyze the matter energy in the model with nonstan- 
dard dynamics. From ([6]), ([8]) and (1201) . we get 



oo 



t\,= I ,h e 2A {<P' 2 + a<p"- l -W 2 ) . (23) 



Finally, using f JT5|) and ( JIB"]) and integrating, we obtain 

E<t> = \ (e 2A(oo) iy(0(oo)) - e 2A (-°°^(0(-oo))) , (24) 

the value of the matter energy for all a. Note that the matter energy depends on the 
asymptotic behavior of the warp factor. 



Now, we follow the same procedure of Ref. [33] and let us focus our study in the case of 
a very small. Thus, the solution of f|T7|) . up to first-order in a becomes 

<P' = \W,~IWI (25) 

and substituting (1251) into (1201) we obtain the potential 



V{<P) = \wl-^Wt-\w\ (26) 



At this point we turn to examine the energy functional (122]) . Substituting the potential 
fl26D in 022]), we get 

E [A, <j>] = |~ dy ll(f~ \W* + 2 - 3 (A> + \W 

+ 1 L dy ( 2 ^ 4 + l w * - + 1 L d 4 ^ • ^ 

From (1271) we can conclude that the solutions of the first order differential equations, ( 1T61) 
and (125]) . are those that minimize the energy functional. Thus, the first order differential 
equations can be seen as the BPS equations and the energy functional could play the role of 
a BPS energy in such scenario. The value of the energy system for a very small is given by 

E [A, 0] = |e 4A(oo) ^(0(oo)) - e 4A( - oo) ^(0(-cx))) I . (28) 



Again, the asymptotic behavior of the warp factor plays a leading role in the value of the 
energy functional. 

Now, we find explicitly the solutions for (I16p and (]2"5]) . which minimizes the energy 
functional. Then, the solution for (1251) becomes 



a 



<Ky) = My) - j w^MvWiMv)) , 

where 4>o(y) is the solution when a = 0. From (Tl6|) and ([29]), we obtain 

A(y)=A (y) + ^W(My)) 2 , 



(29) 



(30) 



>->^" 1-WI--W 2 - — wf 
4 + 3 16 



where Ao(y) represents when a = 0. The matter energy density given by 

p = e 2 

substituting (129]) and (130]) in (13TT) . we obtain 

P = Po - ^ e 2Ao ^ (eW^WjW - 10W 2 Wj + 3W} + ^W 4 



is 



(31) 



where 



fin;- -lu 2 



(32) 



(33) 



4>=4>o 



note that the energy density ( 1321) is a little bit different from that given in Ref. 33 ]. 

To show the validity of the solutions, (129]) . (130|) and ( 132|) . we consider the superpotencial 
W(4>) of the form jg] 

(34) 



where 6 and c are positive parameters. The classical solutions for (129]) and ( 130]) are given by 
<p(y) = \/tt arcsin [tanh (c y)] — ^ v ^ a 5 3 / 2 c 2 tanh (c y) sech (c y) , (35) 



and 



3cv 

A(y) = b In [sech (c y)] + — b 2 c 2 tanh 2 (c y) . 



(36) 



The profiles of the matter energy density is shown in Fig. (pQ) for some values of a. As in Ref. 
331 ] the numerical study gives full support to the analytical expressions for a very small. 
The Fig. ([1]) clearly shows that the brane is localized at y = 0, because this region has a 
positive matter energy density. The contribution of the nonstandard kinetic term modifies 



the profile without altering the symmetrical form of matter energy density, as noted in Fig. 
([1]). It is known that the a-parameter, which indicates the strength of fourth-order kinetic 
term, is used to modify the dynamics of the scalar field contributes to thicker the brane, as 
reported in [33!]. The profiles of the matter energy density and the Ricci scalar are shown 
in Fig. 02]) for a = 0.1. Note that the presence of regions with positive Ricci scalar is 
connected with the localization of the brane and it reinforces the conclusion of the analyzes 
from the matter energy density. Also note that far from the brane, R tends to a negative 
constant, characterizing the AdS§ limit from the bulk, as reported in [27J. A similar behavior 
is obtained for a = 1 and a = 10. The profiles of the warp factor is shown in Fig. ([3]) for 
some values of a. Fig. (|3]) shows that e 2A — > as y — > ±oo, therefore the matter energy 
and the energy of system (128"j) both are zero. 





0.8 






0.6 

,) 






7 


\- 
i 




J 

f().4 

\ 


~\\ 




/: 2 

) ■ . 




^s.^' '' -°- 2 





FIG. 1: The profiles of the energy density for b = 2/3, c = 1, a = 0.1 (dotted line), a = 1 (dashed 
line) and a = 10 (thin line). 




FIG. 2: The profiles of the matter energy density (thin line) and Ricci scalar (thick line) for b = 2/3, 
c = 1 and a = 0.1. 
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FIG. 3: The profiles of the warp factor for b = 2/3, c = 1, a = 0.1 (dotted line), a = 1 (dashed 
line) and a = 10 (thin line). 



Now on, we can analyze the exact solutions in the thin brane limit (c — > oo), we obtain 



v 66 

<P(y) = ^j-Trsgn(y), (37) 

A(y) = -bc\y\, (38) 

from this we can conclude that the braneworld model with nonstandard dynamics treated 
here is consistent, because it reduces to the Randall-Sundrum model (thin brane model). 

III. FERMION LOCALIZATION 

The action for a Dirac spinor field coupled with the scalar fields by a general Yukawa 
coupling is 

S = J (fx^\ ^r M V M * - r]^F((/))^] , (39) 

where 77 is the positive coupling constant between fermions and the scalar field. Moreover, 
we are considering the covariant derivative Vm = du + \ ^tf^A^B > where A and B, denote 
the local Lorentz indices and u^ 3 is the spin connection. Here we consider the field as a 
background field. The equation of motion is obtained as 

iY M V M ^ -r]F{4))^ = 0. (40) 

At this stage, it is useful to consider the fermionic current. The conservation law for J M 
follows from the standard procedure and it becomes 

V M J M = ^ (V M r M ) (41) 



where J M = ^T M <H. Thus, if 

V M r M = 0, (42) 

then four-current will be conserved. The condition (142!) is the purely geometrical assertion 
that the curved-space gamma matrices are covariantly constant. 

Using the same line element ([2]) and the representation for gamma matrices T M = 
(e _j4 7 /i , — 27 s ), the condition ( 1421) is trivially satisfied and therefore the current is conserved. 
The equation of motion (I4"0"j) becomes 

[i-fdp + 7 5 e A («9 y + 2d y A) - 7] e A F(4>)] V = 0. (43) 

Now, we use the general chiral decomposition 

*(x,y) = ^2^L n (x)a Ln {y) + y ^^ Rn {x)a Rn {y), (44) 

n n 

with ipL n (x) = —l br >pL n {.x) and ipR n (x) = rfipR^x). With this decomposition ipL n (x) and 
ijjR n (x) are the left-handed and right-handed components of the four-dimensional spinor 
field, respectively. After applying ( 1441) in ( T43l) . and demanding that i^d^ipLn = f^n^Rn and 
i^d^ipRn = m n ipL n , we obtain two equations for ai n and 

[d v + 2d y A + riF((f>)] a Ln = m n e- A a Rn , (45) 

[d y + 2d y A - r)F{<t>)} a Rn = -m n e- A a Ln . (46) 
Inserting the general chiral decomposition ( 1441) into the action (1391) . using (1451) and (14*61) and 



also requiring that the result take the form of the standard four- dimensional action for the 
massive chiral fermions 

S = J2 I ^ X - m n) ^ (47) 

where ip n = ipL n + ipR n and m n > 0, the functions a^ n and must obey the following 
orthonormality conditions 

/CO 
dye 3A a Lm a Rn = 5 LR 5 mn . (48) 
-00 

Implementing the change of variables 

z= f y e- A ^dy', (49) 
«l„ = e~ 2A L n and a/j n = e~ 2A R n , we get 

+ 14(4^ = 7^, (50) 
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- + V L {z)R n = m 2 n R n 



(51) 



where 



V L (z) = r } 2 e 2A F 2 {^)-r ] d z {e A F{ ( j>)) 
V R (z) = r ] 2 e 2A F 2 ( ( p)+r ] d z (e A F( 



(52) 
(53) 



Using the expressions d z A = e A ( y J_dy A and d z F = e A ^d y F, we can recast the potentials ( 152 
and ( )53l) as a function of y {28]. 29] 



V^fo)) = V L (z{y))\ v ^ v . 



(54) 
(55) 



It is worthwhile to note that we can construct the Schrodinger potentials Vl and Vr from 
eqs. (gU) and fl55|) . 

At this stage, it is instructive to state that with the change of variable (149]) we get a 
geometry to be conformally flat 



ds 2 = e 2A ^ (rj^dx^dx" - dz 2 ) . 



(56) 

Now we focus attention on the condition (j42p for the line element fl56l) . In this case we 
obtain 

V M r M = z(«9 2 ^))e-^) 7 5 . (57) 



Therefore, the current is no longer conserved for the line element f[56|) . It is known that, 
in general, the reformulation of the theory in a new conformal frame leads to a different, 
physically inequivalent theory. This issue has already a precedent in cosmological models 
42]. Recently, other inconsistencies are stated in [30] - 

Under this arguments, we only use the change of variable (1491) to have a qualitative analysis 
of the potential profiles f l54"|) and fl55|) . which is a fundamental ingredient for the fermion 
localization on the brane. 

Now we focus attention on the calculation of the zero mode. Substituting m n = in (j4"5]) 
and (1461) and using ql„ = e~ 2A L n and aR n = e~ 2A R n , respectively, we get 



L oc exp 
Ro oc exp 



-r, / dy'F{(j>) 



V I dy'F{4>) 




(58) 
(59) 
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This fact is the same to the case of two-dimensional Dirac equation 43]. At this point 
is worthwhile to mention that the normalization of the zero mode and the existence of a 
minimum for the effective potential at the localization on the brane are essential conditions 



for the problem of fermion localization on the brane. This fact was already reported in [29]. 

In order to guarantee the normalization condition (I48p for the left-handed fermion zero 
mode (1581) . the integral must be convergent, i.e 



dyexp 



-A(y)-2 V I dy'F(<j>(y')) 
o 



< oo. 



(60) 



This result clearly shows that the normalization of the zero mode is decided by the asymp- 
totic behavior of F((j)(y)). Furthermore, from (l54"j) and (l55j) . it can be observed that the 
effective potential profile depends on the F(cf)(y)) choice. This fact implies that the exis- 
tence of a minimum for the effective potential Vi(z(y)) or Vn(z(y)) at the localization on the 
brane is decided by F(<p(y)). This point will be more clear when it is considered a specific 
Yukawa coupling. Therefore, the behavior of F(<p(y)) plays a leading role for the fermion 



localization on the brane 



29j. Having set up the two essential conditions for the problem 



of fermion localization on the brane, we are now in a position to choice some specific forms 
for Yukawa couplings. 

A. Zero mode and fermion localization 



From now on, we mainly consider the simplest case F(<j)) = <fi. First, we consider the 
normalizable problem of the solution. In this case, from Eqs. (135|) and (1361) the integrand 
in ( 160]) can be expressed as 



/ = exp 



—b In (sech (c y)) — — b 2 c 2 tanh 2 (c y) — r]V6b I{y) — r\ — - — b 3 ^ 2 c sech(cy) 



(61) 

where I = J dy' arcsin [tanh(q/)]. Following the same procedure of Ref. 18], we only need 
to consider the asymptotic behavior of the integrand. It becomes 



/ — > exp 




v K \i ~ — " c I \ y\ 



(62) 



This result clearly shows that the zero mode of the left-handed fermions is normalized only 
for r] > -Wj. Note that the asymptotic behavior of the normalization condition for this 
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case is independent of a. Now, under the change r\ — > —77 (L —> Rq) we obtain that the 
right-handed fermions can not be a normalizable zero mode. The shape of the potentials 
for this case are shown in Fig. (J3J) for some values of a. The Fig. (Jlk) shows that the 
potential of left-handed fermions, Vl, is indeed a volcano- like potential. The shapes of the 
matter energy density, Vl potential and zero mode for this case are shown in Fig. §5§ for 
a = 0.1. A similar behavior is obtained for a = 1 and a = 10. The Fig. ([5]) clearly shows 
that the effective potential Vl has a minimum at the localization of the brane, therefore, 
this clearly shows that the zero mode of the left-handed fermions is localized on the brane. 
On the other hand, the figure (Jlh) shows that the depth of the well structure decreases 
as a increases. From this can be conclude that the ability to trap fermions of the effective 
potential Vl is inversely proportional to a. Figure (jib) shows that the potential Vr is always 
positive. This effective potential has a maximum that decreases as a increases. Therefore, 
none bound fermions with right chirality can not be trapped by the potential Vr. The 
analytic expressions for this model are only valid for a small. For a, which is not necessarily 
small, the numerical study is used. The numerical study done for a large range of values of 
a bear out our results. 




FIG. 4: Potential profile: (a) (V l (v))a (left) and (b) (V R (y)) A (right) for 77 = 1, b = 2/3, c = 1, 
a = 0.1 (dotted line), a = 1 (dashed line) and a = 10 (thin line) 



IV. CONCLUSIONS 

We have reinvestigated the braneworld model constructed out from one scalar field with 
nonstandard kinetic terms coupled with gravity. We have considered the model C = K(X) — 



V((f)), where K = X + a\ X\X (type I model in Ref. 33] ). This work completes the analysis 
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FIG. 5: The profiles of the energy density (thin line), (Vl)a (thick line) and zero mode (dashed 
thick line) for b = 2/3, c = 1, rj = 1 and a = 0.1. 



of the research in Ref . 33| . We showed that the equations of motion for all a can be deduced 



from the functional E[A, <ft] (1221) . as done by Townsend {4] in the case of standard dynamics. 
Furthermore, we showed that for a small the solutions of the first order differential equations, 
( I16p and ( )25l) . are those that minimize the energy of system. In contrast, for a, which is 
not necessarily small, the solutions of the first order differential equations not necessarily 
minimize the energy of system. Also, we showed that the value of the matter energy and 
the energy of system depends on the asymptotic behavior of the warp factor. We found an 
expression for the matter energy density that differs slightly from [33]. The numerical study 
gives full support to our matter energy density expression for a small. 

We also have investigated the localization problem of fermions for the type I model. We 
have used the simplest Yukawa coupling r/\i/0\l/ between the scalar and the spinor fields. 
In order to guarantee the normalization condition for the zero mode, we showed that the 
zero mode of left-handed fermions is normalizable under the condition ri > - a/tt and it is 
independent of a. For this kind of solution, the effective potential of left-handed fermions 
Vl is a volcano-like potential. Vj, has a minimum at the localization of the brane (y = 0), 
therefore the zero mode of the left-handed is localized on the brane. On the other hand, 
the value of a adjust the minimum of Vl, the depth of the well structure decreases as a 
increases. Therefore, we can conclude that the ability to trap fermions of Vl is inversely 
proportional to a. The right-handed fermions can not be localized on the brane, this fact is 
a consequence of the absence of a normalizable zero mode. For a not necessarily small, the 
numerical study done for a large range of values of a bear out our results and conclusions. 
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Additionally, we showed that the change of variable dz = e~ A ^dy leads to a non con- 
served current, because the curved-space gamma matrices are not covariantly constant. The 
effects of non-conserved current in the issue of massive modes is currently under considera- 
tion and will be the subject of another thorough study. 
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